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Abstract 

We extend the discrete Regge action of causal dynamical triangulations to include discrete versions of 
the curvature squared terms appearing in the continuum action of (2 + l)-dimensional projectable Hofava- 
Lifshitz gravity. Focusing on an ensemble of spacetimes whose spacelike hypersurfaces are 2-spheres, we 
employ Markov chain Monte Carlo simulations to study the path integral defined by this extended discrete 
action. Wc demonstrate the existence of known and novel macroscopic phases of spacctimc geometry, 
and we present preliminary evidence for the consistency of these phases with solutions to the equations 
of motion of classical Hofava-Lifshitz gravity. Apparently, the phase diagram contains a phase transition 
between a time-dependent de Sitter-like phase and a time-independent phase. We speculate that this 
phase transition may be understood in terms of deconfinement of the global gravitational Hamiltonian 
integrated over a spatial 2-sphere. 
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1 Connections 



An intriguing body of evidence hinting at a deep connection between Hofava-Lifshitz gravity^ and causal 
dynamical triangulations has recently accumulated in the literature. First, one of us demonstrated consis- 
tency of the spectral dimension computed in Hofava-Lifshitz gravity with the spectral dimension measured in 
causal dynamical triangulations in 3 -I- 1 dimensions |32| . Benedetti et al then verified within causal dynam- 
ical triangulations |19| the prediction of Hofava-Lifshitz gravity for the behavior of the spectral dimension 
in 2 -|- 1 dimensions [32j . Next, Ambj0rn et al exhibited the remarkable resemblance of the phase diagram 
of causal dynamical triangulations [2] to the phase diagram of Lifshitz matter systems |29] as exemplified by 
the Lifshitz scalar field f34|. These authors and one of us also conjectured [21 [34] that the apparent tricritical 
point of the former phase diagram could correspond to the tricritical limit of Hofava-Lifshitz gravity with 
dynamical critical exponent z equal to the dimension of space. Then, both Benedetti et al and Ambj0rn et 
al noted the compatibility of certain solutions of Hofava-Lifshitz gravity with the minisuperspace model fit 
to the expectation value of the geometries emerging from causal dynamical triangulations [2 119) . Recently, 
Sotiriou et al successfully fit the behavior of the spectral dimension of causal dynamical triangulations at 
intermediate scales to a dispersion relation derived from Hofava-Lifshitz gravity [48j . Just days ago, Budd 
argued that the kinetic term in the semiclassical effective action of causal dynamical triangulations exhibits 
a Hofava-Lifshitz-like form [21]. 

This mounting evidence motivated us to extend the Regge action — or, more precisely, the discrete path in- 
tegral measure — used in causal dynamical triangulations to include dependence on the broader class of terms 
appearing in the action of Hofava-Lifshitz gravity. Our purpose is threefold: firstly, to test the applicability 
of causal dynamical triangulations to modified classical theories of gravitation; secondly, to explore quan- 
tum Hofava-Lifshitz gravity with nonperturbative techniques; and, thirdly, to illuminate further the links 
between Hofava-Lifshitz gravity and causal dynamical triangulations. As an initial step toward these goals, 
we have begun to investigate an appropriate reduction of (2 -|- l)-dimensional projectable Hofava-Lifshitz 
gravity using causal dynamical triangulations. This model, though removed from the phenomenologically 
interesting case of 3 -I- 1 dimensions, serves as a simplified yet nontrivial beginning for our research. 

After briefly introducing the formalisms of Hofava-Lifshitz gravity and causal dynamical triangulations in 
section [2| we explain our adaptation of causal dynamical triangulations to Hofava-Lifshitz gravity in section 
[3] We present the results of our initial numerical studies — evidence for the existence of extended phases of 
geometry, the structure of these phases, and their relation to the classical solutions — in section|4j Finally, we 
summarize our conclusions and discuss ongoing and future work in section [5j To streamline our presentation, 
we relegate to appendix |X] the construction of the relevant classical solutions of Hofava-Lifshitz gravity and 
to appendix [B] certain geometric properties of causal dynamical triangulations. 



2 Background 

2.1 Hofava-Lifshitz Gravity 

Hofava-Lifshitz gravity is a field theory of the dynamical metric on spacetime manifolds that carry the 
additional structure of a preferred foliation by spacelike hypersurfaces. The relevant set of gauge symmetries 
is then the group DiSpi-M) of diffeomorphisms of the spacetime manifold Ai that preserve the preferred 
foliation F. The foliation structure leads to one important novelty: the possibility of anisotropic scaling 
with a nontrivial dynamical scaling exponent z measuring the degree of anisotropy between space and time. 
In Minkowski spacetime the relativistic scaling is thus replaced by the anisotropic scaling 

t — > i=bH (2.1a) 
X — > X = 6x (2.1b) 

for constant 6 > 0. Starting with z > 1 at short distances markedly improves the ultraviolet behavior of the 
theory, potentially rendering it power-counting renormalizable. 

^We have decided, by a majority vote among the coauthors of this paper, to foUow the terminology commonly accepted in 
the literature and refer to gravity models with anisotropic scaling as Hoi'ava-Lifshitz gravity, despite the dissenting vote of one 
of us (PH). 
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One might be tempted to rewrite this theory in a manifestly relativistic fashion by integrating in the 
remaining part of the group Diff(A^) of fuh spacetime diffeomorphisms. This procedure leads to a reformula- 
tion of the model as a specific scalar-tensor theory with higher-derivative interactions in which the spacelike 
hypersurfaces of constant scalar field dynamically determine the leaves of the foliation F. Such a relativistic 
rewriting is only equivalent to the original nonrelativistic formulation at the classical level with, moreover, 
subtle regularity conditions on the scalar field's dynamics. At the quantum level the relativistic rewrit- 
ing immediately besets the theory with the notorious problem of time, whereas the nonrelativistic nature 
of the original formulation of Hofava-Lifshitz gravity potentially renders the spacetime metric's dynamics 
more directly compatible with quantum mechanics. Additionally, as we shall show, the original formulation 
of Hofava-Lifshitz gravity is nicely suited to the framework of causal dynamical triangulations, which also 
utilizes a preferred foliation structure in the microscopic definition of the path integral for gravity. 

In a local smooth coordinate chart (t,x) adapted to the preferred foliation F, 'DiSp{J\4) consists of all 
reparametrizations of the form 

t — > f{t) (2.2a) 
X — > x = C{t,x) (2.2b) 

for arbitrary functions / and C. Note that, although reparametrizations of the space coordinates may be 
time-dependent, reparametrizations of the time coordinate must be space-independent. Given its preferred 
foliation, the structure of Hofava-Lifshitz gravity is naturally discussed in the Arnowitt-Deser-Misner for- 
malism fl8L In this formalism the spacetime metric tensor g(t,x) is decomposed into the metric tensor 
7(i,x) on a spacelike hypersurface E of constant time coordinate t, the shift vector N(i,x), and the lapse 
function N{t, x) such that one can reassemble the standard line element as 

ds^ = -iV2(i,x)dt2 +7,j (t,x) [dx' + N\t,yi)dt] [dx^ + {t,:>c)dt] . (2.3) 

Note, however, that recombining the spatial metric tensor, the shift vector, and the lapse function into 



the spacetime line element (2.3) is somewhat misleading since, in the regime with z > 1, different terms 
contributing to ds^ carry different scaling dimensions. 

The shift vector and the lapse function play the role of gauge fields associated with DiSp{A4). Since 



the time coordinate reparametrizations (2.2a) are independent of the space coordinates, a natural choice 
is to restrict the corresponding gauge field N{t,x) also to be only a function of the time coordinate. We 
choose to make this restriction, yielding the so-called projectable version of the theory. (See, for instance, 
[311 1201 WE[ 1521 153) . ) Of course, we have interest in studying the more general nonprojectable version in which 
the lapse function is a spacetime field and of which the projectable version is a dynamical limit. The gauge 
symmetries then permit new terms in the action involving spatial derivatives of the lapse function; finding 
the appropriate realization of such terms in the framework of causal dynamical triangulations is beyond the 
present work's scope. We instead focus on the projectable version of Hofava-Lifshitz gravity, which has a 
particularly clear translation into the language of causal dynamical triangulations. 

Employing the Arnowitt-Deser-Misner decomposition of the metric tensor, we now construct the action 
of projectable Hofava-Lifshitz gravity. We aim to build a power-counting renormalizable unitary classical 
theory of gravitation. Even in relativistic gravity, including higher curvature terms can render the theory 
renormalizable |261 150] : typically, however, such terms come at the cost of sacrificing perturbative unitarity 
and propagating unphysical degrees of freedom |49| . These issues stem from the inclusion of higher temporal 
derivative terms, which necessarily accompany higher spatial derivative terms in a relativistic theory. By 
only permitting higher spatial derivatives, we can in principle avoid problems with perturbative unitarity 
and unphysical degrees of freedom. Working in the Arnowitt-Deser-Misner formalism allows for the rather 
straightforward inclusion of higher spatial derivatives and exclusion of higher temporal derivatives as we 
desire in constructing the action of Hofava-Lifshitz gravity. 

We build the general action in d spatial dimensions as the sum of a kinetic term quadratic in temporal 
derivatives and a potential term of mass dimension 2d in spatial derivatives. The dynamical critical exponent 
z thus equals d. Accordingly, we form the kinetic term from invariants of dtj{t,x.). This derivative alone 
is not covariant under DiSpiM), but the extrinsic curvature tensor K(t,x) of a spacelike hypersurface S 
satisfies this criterion. With K(t,x) having components 

K,j{t,x) - [^tl^J{t,^) ~ V,7V,(t,x) - V,iV,(t,x)] (2.4) 
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for the covariant derivative V associated with "f{t, x), the most general kinetic term invariant under T)iSp(Jli) 
is 

/ dtd''x./W^)N{t) [K,,{t,^)K^^it,^) - XK^t,^)] . (2.5) 

lOTTLz Jj^ 

Here the parameter A arises from the generahzed DeWitt supermetric compatible with the theory's gauge 
symmetries, and K{t,x) is the trace of the extrinsic curvature tensor K(t,x) [501 131j . 

We form the potential term from invariants of 7(t,x) and its spatial derivatives. The most general 
potential term is 

-J— [ dtd''x,/W^)Nit)VMt,^)], (2.6) 
where V Yf{t, x)] is a scalar functional of the spatial metric tensor and its spatial derivatives up to order 2d. 



Putting together (2.5) and (2.6), the action of projectable Hofava-Lifshitz gravity is 



■>HL\ 



g(<,x)]-^ / dtd'^a;v/7(i,x)7V(t) {X,,(t, x)i^'^(i, x) - XK^t,^) - V^[7(i,x)]} . (2.7) 



M 



Note that the coupling constant G is related but not equal to the Newton constant Gm'- the low energy 
Newton constant Gjv is determined by rescaling the time coordinate by an effective speed of light factor 
selected so that the term linear in the spatial Ricci scalar within y[7(t,x)] is correctly normalized with 
respect to the kinetic term [301 [SI]. This rescaling also dictates how the constant term in V [7(i,x)] relates 
to the low energy cosmological constant A. 

For d> 2, with the anticipated z = d scaling at short distances, there is a proliferation of marginal and 
relevant contributions to the potential term; for d = 2, however, the most general action containing only the 
marginal and relevant terms for z = 2 is quite compact: 

SHL[g{t,^)] = I dtd^x^-i{t,^)N{t) [K,,{t,^)K'\t,^) - XK^t,^) 

-aRl{t, x) + PR2{t, x) - 2A] (2.8) 

with i?2(i,x) the Ricci scalar of the spatial metric tensor 7(i,x) and the coupling constant A related to the 
low energy cosmological constant A by the rescaling described above. The potential term's relative simplicity 
in (2 + 1) -dimensional projectable Hofava-Lifshitz gravity further motivates our initially studying this case. 



The equations of motion (A.2) derived from variation of the action (2.8) with respect to the spatial 
metric tensor are independent of the coupling constant /3: for d = 2 the i?2 term is a total derivative, and 
the Gauss-Bonnet theorem determines its integral over a spacelike hypersurface S of the preferred foliation 



in terms of the Eulcr number of S. Similarly, the equations of motion ( A.l ) obtained from variation of the 



action (2.8) with respect to the shift vector are also insensitive to /3. The coupling constant fi does, however. 



appear in the equation of motion derived from variation of the action (2.8) with respect to the lapse function. 
Since the lapse function only depends on the time coordinate, its equation of motion takes the form of a 
spatially integrated Hamiltonian constraint 

•H_L = (2.9) 

for the global Hamiltonian 



Hi = 



/ d^xy/"i{t,^) [K,j{t,x.)K'i{t,x) - XK^{t,x) + aRl{t,yi) - /3i?2(i,x) + 2A] . (2.10) 



At this stage there is an additional choice to make beyond that of restricting to projectable Hofava- 
Lifshitz gravity. Recall the distinction that arises between noncompact and compact spatial topology when 
imposing the local Hamiltonian constraint in relativistic theories of gravity. In the former case, after imposing 
appropriate asymptotic fall-off conditions and accounting for possible boundary contributions, the zero mode 
'H_L of the local Hamiltonian constraint does not vanish on physical states. Instead, this mode is one of the 
conserved asymptotic charges, namely the total energy. In the latter case the zero mode T-L± of the local 
Hamiltonian constraint vanishes on all physical states since the relativistic gauge symmetries mix it together 
with the other gauge symmetry generators. 
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Figure 1: From left to right: the (3, 1)-, (2,2)-, and (1, 3)-tetrahedra. The abscissa indicates the number of 
vertices on an initial triangulated spacelike hypersurface, and the ordinate indicates the number of vertices 
on the next triangulated spacelike hypersurface. 



Now, in Hofava-Lifshitz gravity, as pointed out in [33J [Ml |S3], this situation is somewhat different: 
while the case of noncompact spatial topology remains essentially unchanged, a novelty arises in the case 
of compact spatial topology. Imposing the vanishing of Ji^ on physical states is not mandatory: either 
one may treat as a gauge symmetry generator and thus impose the condition "Hj^ = 0, or one may 
treat H± as a conserved global charge measuring the energy levels of physical states in the Hilbert space. 
Of course, the second option requires us to check that this Hamiltonian's spectrum is bounded from below 
on physical states. This novel situation was explicitly encountered in [54] whose authors constructed from 
bosonic systems on a rigid lattice the renormalization group fixed points corresponding to the free field limit 
of Hofava-Lifshitz gravity with z = 2 and 2 = 3. There the specific choices that one could make at the level 
of the microscopic lattice Hamiltonian gave rise to the option of imposing or not imposing the condition 
H± = for the effective low energy gravitons. 

) ) in section [3] for use as a functional on simplicial manifolds, we find 
as a conserved global charge — appears more natural in the setting of 



When we discretize the action (2. 
that the second option 



-treating H± 

causal dynamical triangulations. The geometric restrictions derived by employing simplices of fixed edge 
lengths suggest that the lapse function has been effectively set to a constant, though without leading in any 
discernible way to the imposition of the constraint (2.9 1. Also, the i?2 term with coupling constant jS in 
the action (2.8) does not affect any observables for the topological reasons stated above. The value of H± 



does, however, depend on (3, which suggests that enforcing the constraint (2.91 is not completely consistent 



In section 4.3 where we compare the geometries emerging from our numerical simulations to the classical 



solutions for the constraint (2.9) not imposed, we find evidence supporting our selection of the second option. 



These comparisons are, moreover, not without weight: as briefiy discussed in j34j and further demonstrated 
in [33], the classical phase diagram for the theory with variable lapse differs significantly from that of the 
present setting. The classical theory that we quantize in section |4] is thus defined by the action (2.81 for 
fixed lapse. 



2.2 Causal Dynamical Triangulations 

The causal dynamical triangulations approach aspires to define the continuum limit of a quantum theory of 
gravity by appealing solely to those nonperturbative tools applied with great success to the quantization of 
local gauge field theories. (See 38, 39 for reviews.) In particular, the approach invokes a lattice regularization 
of the path integral for gravity and then utilizes the renormalization group to explore its continuum limit. 
In causal dynamical triangulations one thus attempts to define this path integral as the partition function of 
a statistical model of dynamical geometry. (See |:3i. il6] for recent reviews.) As in several previous programs 
of lattice quantization of gravity, the statistical ensemble of geometries is comprised of simplicial manifolds 
triangulated by a fixed set of {d + l)-simplices. In 2 -I- 1 dimensions this set consists of those 3-simplices or 
tetrahedra pictured in figure 1. 

The novel feature of causal dynamical triangulations lies in the imposition of an additional restriction on 
the geometries permitted to enter the path integral: these simplicial manifolds must possess a global foliation 
by spacelike hypersurfaces of constant discrete time. The global foliation was originally introduced to enable 
a Wick rotation from Lorentzian to Riemannian signature. The causal structure of a Lorentzian triangulation 
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could thus be faithfully retained in its Wick rotated version, only the latter triangulation being amenable to 
numerical analysis [71IH1E]. A simplicial manifold is endowed with the foliation structure as follows: every 
spacelike hypersurface, all of a chosen fixed topology, is triangulated with equilateral d-simplices, and then 
the vertices of adjacent spacelike hypersurfaces are connected by timelike edges so as to produce only the 
{d + l)-simplices of the fixed set. Spacelike edges have length squared Ig^ = c? , defining a lattice spacing 
for the triangulation, and timelike edges have length squared = —'qa? for 77 > [71 |H] . Note that our 
parameter t] corresponds to the parameter a typically used in the causal dynamical triangulations literature. 
More concretely, one aims to approximate the path integral for general relativity. 



Z[7(ty,x)|7(i„x)]= / " I?g(i,x)e*^--fe(*''')l, (2.11) 

■^7(t.,x) 

by the path sum 

^[r/|r.] = E7^^"^"'^- (2-12) 
r 

In (2.11) the path integral is taken over all physically distinct spacetime metric tensors g(t,x) interpolating 
between the initial and final boundary geometries specified by the spatial metric tensors 7(i,;, x) and 7(t/, x), 
and Seh denotes the Einstein-Hilbert action 

5Eff[g(t,x)] = — ^ / AtA'^x ^-g{t,^) [i?(t,x)-2A]. (2.13) 



In (2.121 the path sum is taken over all causal triangulations T interpolating between the initial and final 



boundary geometries specified by the triangulations and T f [H [8] , and Sr denotes the Regge action 



SttGn 16ttGn ^ 

her seT 



Here, /i is a 1 — 2)-dimensional hinge having area and deficit angle 5h, and Vs is the spacetime volume 
of a (d + l)-simplex |51) . The measure factor ^ is the inverse of the order of the automorphism group of 
7", included to account for discrete symmetries arising in the structure of T. 

Before proceeding further with our discussion of causal dynamical triangulations, we make a few clarifying 



remarks on the implications of the above prescription for regularizing the continuum path integral (2.11 1. The 
restriction to simplicial manifolds admitting a preferred foliation effectively changes the path integral measure 
and integration domain. Even though we continue to employ the Regge action, these changes may result 
in our model belonging to a universality class different from that of Euclidean dynamical triangulations. In 



1 + 1 dimensions, where the path sum (2.12 ) can be evaluated analytically, this is known to be the case [T]: by 
preventing the birth of baby universes, the foliation requirement alters the critical exponents, indicating that 
Euclidean and causal dynamical triangulations occupy different universality classes. In higher dimensions 
the same situation seems very likely to hold true; otherwise, we should expect our model to fall into one of 
the universality classes in which no smooth macroscopic limit exists, and the spacetime geometry exhibits a 
branched polymer or crumpled behavior. 

The change in the path integral measure and integration domain may also translate into a change in our 
theory's effective action such that it no longer assumes the general relativistic form. Recall that, although 
we may locally treat the choice of a preferred foliation as a gauge choice, we cannot generally make such 
a choice globally. A gauge choice typically introduces a Jacobian into the path integral measure, and a 
mismatch between the path integral measure and the gauge choice can lead to a breaking of the gauge 
symmetry. Starting with the Regge action, the question of whether counterterms sensitive to the preferred 
foliation are generated thus remains open. If such counterterms appear, then the resulting effective action 
likely corresponds to Hofava-Lifshitz gravity at some particular values of its couplings. The Regge action's 
bare coupling constants Gjv and A are only indirectly related to the continuum renormalized values of the 
Newton constant and the cosmological constant at long distances, so we must not naively identify the former 
with the latter. 



Continuing with our discussion of causal dynamical triangulations, the action (2.141 simplifies consider- 



ably for the simplicial manifolds contributing to our ensemble since only a small set of simplices with fixed 
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geometries is used to construct them. Still, the path sum (2.12 ) has resisted all attempts at analytical compu- 



tation for d > 1, so exploration of its properties has primarily employed numerical techniques. In particular, 
Markov chain Monte Carlo methods are used to simulate the path sum (|2.12|) 4, 5, 9, 10, llt lT2l[T4l [T5 , .24 , .40). 



To render the path sum (2.121 amenable to such analysis, we must first Wick rotate the real time action 
( 2.14[ ) to imaginary time, a well-defined process owing to the foliated structure of causal dynamical triangu- 
lations. We perform the Wick rotation by analytically continuing the parameter ry in the lower half complex 
plane [Hj . For our case of interest — 2 + 1 dimensions with spacelike hypersurfaces having the topology of 
and periodic time coordinate having the topology of the 1-sphere — the Regge action becomes 



(2.15) 



after Wick rotation and application of various topological relations for 77 = 1 [3]. Here, Nq is the number 
of vertices and N3 is the number 3-simplices in the triangulation T. The coupling constants fco and are 
related to the bare couplings Gn and A as 



fco 



4Gn 

a^A a 



We have thus transformed the path sum (2.12) into the statistical partition function 



T 



1 



(2.16a) 
(2.16b) 

(2.17) 



Computer simulations of the partition function (2.17) in both 2 + 1 and 3 + 1 dimensions have thus 



far provided considerable support for the existence of an extended phase of geometry possessing not only a 
semiclassical limit on large scales, but also a quantum regime on small scales [H [HI UHl [HI EH ESI Ell EH EH 
mi 130] . In particular, the average observed geometry matches well that of (possibly deformed) Euclidean de 
Sitter spacetime at both the classical and semiclassical levels [lliniEniEIlE2E31EllII3Enil21]- Furthermore, 
studies of the spectral dimension of this extended phase of geometry have revealed an apparent dimensional 
reduction to effective 2-dimensionality on small scales dSlEHlSD]- This phenomena of dynamical dimensional 
reduction — particularly, the extrapolated value of the minimal dimensionality — has elicited comparisons of 
causal dynamical triangulations to both the asymptotic safety approach and Hofava-Lifshitz gravity. In 
the former theory the effective change in the spectral dimension apparently results from large anomalous 
dimensions near the nontrivial fixed point, while in the latter theory the spectral dimension flows to the 
short distance value of 2 as a result of the anisotropic scaling near the Gaussian fixed point [32l |4^ . Further 
evidence for such dimensional reduction has also surfaced in other approaches to the quantization of general 
relativity |23j. 

As in Euclidean dynamical triangulations, additional phases of geometry also emerge |[43 . At first 
these phases were viewed as unphysical, but current interpretations favor their role as further phases in the 
vicinity of a multicritical fixed point Specifically, in both 2 + 1 and 3 + 1 dimensions there exists a 
phase characterized by spacelike hypersurfaces that effectively decouple from one another [H EH SD] > and in 
3 + 1 dimensions a second additional phase, distinguished by its large Hausdorff dimension, appears [THHO] . 
Recently, Ambj0rn et al have argued that the transition between this last phase and the physical phase 
is of second order 0. This finding raises the possibility of rigorously defining a continuum limit of causal 
dynamical triangulations. 



3 A Discrete Action for Horava-Lifshitz Gravity 



We now derive a discrete form of the action (2.8 ) suitable for analysis with the techniques of causal dynamical 
triangulations. As above we assume a topology of x S^, primarily motivated by the relative ease of 
numerically analyzing such compact spacetimes. As in the lattice regularizations of local quantum field 
theories [381 139j. the precise details of the discretization do not matter since universality ensures that many 



of the details at the lattice spacing scale become irrelevant in the long distance limit. Our primary goal in 
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constructing a discrete analogue of the action (2.8 1 is thus to build an action sufficiently specific so that the 
continuum limit lies in the universality class of Hofava-Lifshitz gravity yet sufficiently generic so that the 
continuum limit does not depend on all of the renormalizcd coupling constants in this universality class. Of 
course, we can only judge whether or not we have achieved these goals once we have thoroughly studied the 
quantum theory of the model defined below. 



Working from the continuum action (2.8), we allow two technical criteria to guide our construction of 
its discrete version: first, the discrete action should manifestly reduce to the Regge action used in causal 
dynamical triangulations when the bare coupling constants A and a assume their general relativistic values, 
and, second, the transfer matrix corresponding to the discrete action defined on the space of boundary 
geometries should yield a well-defined Hamiltonian. These two criteria apply to the classical discrete action 
constructed below. In the quantum theory defined via the path integral based on this action, we naturally 
anticipate a nontrivial relationship between the bare coupling constants and the renormalized coupling 
constants. Depending on this relationship, the quantum theory's long distance continuum limit may or may 
not be relativistic. 

To implement the first criterion, we use the Gauss-Codazzi equation, 

R{t, x) = i?2 (t, x) - [K^ (t, x) - Kij {t, x)K'^ (t , x)] + Total Derivative, (3.1) 



to rewrite the action (2.8) as 

Shl[sM] = / dtd^x^-g{t,x) [i?(t,x)-2A] 

167rG Jj^ 



1- A 
"l67rG 



M 

dt d^xy/-/{t, x)N{t)K^{t, x) 



M 



IGnG 



dtd^xy/'y{t,x.)N{t)Rl{t,x) (3.2) 

M 



up to irrelevant boundary terms. We neglect the term in the action (2.8) with coupling constant f3 as it 



contributes only an additive constant for the spacetime manifolds under consideration. In this form the 



action (3.2) straightforwardly reduces to the Einstein-Hilbert action (2.13) when the coupling constants A 



and a take on their general relativistic values of one and zero. We use the discrete action (2.15) for the 



Einstein-Hilbert portion of the action (3.2 1, leaving us the task of determining the discrete analogues of the 



and i?2 terms. We construct these terms below, postponing discussion of the second criterion. 
3.1 Curvature Squared Terms in Regge Calculus 

There exist well-established prescriptions for constructing the Ricci scalar and the trace of the extrinsic 
curvature tensor in Regge calculus. The former is defined in terms of deficit angles 6h about (d -I- 1 — 2)- 
dimensional hinges h of a {d+ l)-dimensional simplicial manifold [51] : 

f dtd^x ^-g[t,x)R{t,x) = 2Y,Ah5h- (3.3) 
•'M her 



The latter is defined in terms of angles Tph between the normal vectors to the two d-simplices intersecting at 
the (d + 1 — 2)-dimensional hinge h within the spacelike hypersurface E of a (d -f l)-dimensional simplicial 
manifold 

j d'^x ^-i{t,x)K{t, x) = ^ AhiPh- (3.4) 
•'^ he's 
Technically, these curvatures have a distributional definition on the hinges, which complicates the construc- 



tion of curvature squared terms. In particular, simply taking the square of (3.3) or ( |3.4[ ) to define the discrete 
versions of the i?2 o'' terms leads to a mathematically ill-defined continuum limit [6]. Accordingly, we 
adhere to the philosophy of [6l I27j when building discrete analogues of curvature squared terms, adopting 
the alternative scheme of volume sharing. As discussed in [B] , we make the identification 

d'xyWx)^^t,x) ^o^'^C-^X (3-5) 
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for a curvature scalar ^(i,x), where o is the object assigned the curvature density, Or{T) is the set of all 
objects o on the spacelike hypersurface E labelled by discrete time coordinate r in the triangulation T, Vo is 
the appropriate volume of the object o, and Vo is the share-volume of the object o, namely the volume of 
all top-dimensional objects containing o. Using this scheme we now address the R2 and terms in turn. 

3.1.1 Rl Term 

For the term there are clear choices for the objects o and the top-dimensional objects containing o: since 
the Ricci scalar characterizes the intrinsic geometry of a 2-dimensional spacelike hypersurface, the objects o 
are vertices v and the top-dimensional objects are spacelike triangles A. This is completely consistent with 
the usual prescription for the Ricci scalar in Regge calculus. For a vertex v the deficit angle is 



Sy ~ '2,7: 



(3.6) 



for the number N^(v) of spacelikc triangles containing v, the volume is 

K = 1, 

and the share- volume is 



V3 2 



(3.7) 



(3.8) 



Add 



since all of the spacelike triangles are equilateral. Employing (3.5) and noting that the most natural dis- 
cretization of the time integral is 

(3.9) 



/rja, 



we make the identification 

dt d^xy/ j{t,x)N{t)Rl{t,x) 



M 



E E 



a Na{v) 



(3.10) 



up to multiplicative factors, where Vr{T) denotes the set of vertices v belonging to the spacelike hypersurface 
labeled by the discrete time coordinate r. 



3.1.2 K'^ Term 

For the term the choices of objects o and top-dimensional objects containing o are not as clear. The 
extrinsic curvature captures in part how the spacelike hypersurface is embedded in the spacetime manifold. 
The discrete analogue of the term must reflect this geometric information, requiring that we appropriately 
account for how tetrahedra connect to the spacelike hypersurface. This observation suggests that we take 
tetrahedra as the top-dimensional objects contributing to the share-volume. Furthermore, in the continuum 
K"^ scales as an inverse length squared, implying that the objects o be spacelike triangles. We thus need to 
assign a deficit angle (5a to a spacelike triangle. We largely follow the treatment of [37]. Consider a spacelike 
hypersurface of the triangulation T. For each spacelike triangle we may define a future-directed normal 
vector at the center of the (3, l)-tetrahedron of which it forms the base. A natural measure of extrinsic 
curvature at the common edge e between two adjacent such spacelike triangles is the deficit of the angle 
that the normal vector traces out as it is parallel transported from its own (3, l)-tetrahedron to the adjacent 
(3, l)-tetrahedron. In figure 2 we illustrate this construction for the particular case of three (2, 2)-tetrahedra 
stacked between the two adjacent spacelike triangles. In general, the deficit angle for the common edge is 

,5, = i(vr-20f^'-0f^)iV(\,)(e)) (3.11) 

with 9)^ ' and 9)^ ' the Lorentzian dihedral angles about spacelike edges and 7V^2 2) (e) the number of future- 
directed (2, 2)-tetrahedra attached to the common edge e. We give the values of 0]^ ' and 9\^' in appendix 
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Figure 2: An embedding in three dimensions of two (3, l)-tetrahedra (solid black) joined by three (2, 2)- 
tetrahedra (thin black) all sharing a common edge. A vector perpendicular to the triangular base of a (3, 1)- 
tetrahedron rotates through an angle tt — 26^^'^'' — 39^^'^^ as it is parallel transported across the common 
edge. 



[B] We thus assign to a spacelike triangle the deficit angle 



3A 



- (Stt 



66 



(3,1) 



(2,2) 



(3.12) 



which, we note, is reminiscent of a trace. The volume of a spacelike triangle is 



(3.13) 



and the share volume is 



4^^3,1) 



+ e^'4,2)(A) 



(3.14) 



with V£ ' ' and ' the Lorentzian 3-volumes of the respective tetrahedra. We give the values of 
and V^'^'^-* in appendixjB] The share volume assumes this value since a given spacelike triangle has four (3, 1)- 
tetrahedra and 7V^^2 2)(^) (2, 2)-tetrahedra in its immediate future. According to the prescription (3.5), we 
make the identification 



Stt- 661 



(3,1) 



4F, 



(3,1) 



+ ^i'''^<2.2)(A) 



(3.15) 



up to multiplicative factors, where Tf'^[T) denotes the set of spacelike triangles A belonging to the spacelike 
hypersurface labelled by the discrete time coordinate r. 

The discretization (3.15) of the term does not, however, respect our second criterion. Following [8], to 
ensure the existence of a well-defined Hamiltonian on the space of boundary geometries, we must make ( 3.15 ) 
time-reversal invariant. A straightforward calculation shows that this invariance guarantees positivity of the 
squared transfer matrix, which, along with the transfer matrix's symmetry, yields a well-defined Hamiltonian 
[5|. To realize time-reversal invariance, we add an analogous term for past-directed (2, 2)-tetrahedra. The 
complete discrete analogue of the term is 



E E «^ 



Stt - 661 



(3,1) 



9(2,2)^T 



(2,2) 



(A)| 



Stt - 66* 



(3,1) 



(2,2) 



(A)| 



4F, 



(3,1) 



^f''^<2,2)(^) 



4y, 



(3,1) 



(3.16) 



Note that we did not require such an adjustment for the i?2 term since it depends only on the intrinsic 
geometry of the spacelike hypersurface. 
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3.2 Imaginary Time Action 



Putting together (3.10) and (3.16), our discrete action for Hoi^ava-Lifshitz gravity becomes 



Shl\T] — Scdt[T] 
1- A 



T AeTSi(r) 



\37r-69 



(3,1) fl(2,2) 



^('2,2) (A) I 



(3,1) 



+ ^i'''^4,2)(A) 



E E 



/v{6-Na{v)Y 



IGttG^ ^ a Na(v) 



377-66* 



(34) ^(2,2) 



%2)(A)| 



(3,1) , -,.(2.2) 



+ ^(2.2)(A) 



(3.17) 



Wick rotating the action (3.17) to imaginary time, we find that 



1 - A 
167rG 



E E «^ 



3.-6er^-e^)4^,)(A))^ (3.- 



(2,2) 



+ 



167rG' 



E E 



«evv(r) 



V^ (6-iVA(f))^ 



66* 



(3.1) _ ^(2 



4y, 



(3,1) 



+ V^i'''^<2)(A) 



(3.18) 



We give the values of ' , 0^ ' ^\ V^''^\ and 14^'^^ in appendix |b] We use the action (13. 181) for a = 1 and 
77 = 1 in the path integral quantization analyzed below. The value of the length a has no a priori physical 
meaning, and any value of the parameter 77 > ^ is permitted. 



4 A Lattice Quantization of Horava-Lifshitz Gravity 
4.1 Numerical Implementation 

Following the analyses of[l|51il[inilIIl[nilIllI3[Mllin],we explore the partition function 

^[r/|r.] = E7^^"""'^' (4-1) 

employing Markov chain Monte Carlo methods. In our simulations so far, we have fixed the topology of 
the spacelike hypersurfaces to be that of and have fixed the total number T of spacelike hypersurfaces, 
introducing a discrete time coordinate r that enumerates these hypersurfaces. In the simulations reported 
below, we have set T = 64. Additionally, we impose periodic boundary conditions on this time coordinate, 
endowing it with the topology of . Consequently, the initial triangulated spacelike hypersurface Vi and the 
final triangulated spacelike hypersurface of each causal triangulations T entering the partition function 



(4.1) are identified. 

Next, we hold the number of tetrahedra in each triangulation T approximately fixed; otherwise, the 



weights e ^hl^'T^ appearing in the partition function (4.1 ) may grow without bound, eventually causing our 



computer code to crash. To implement the constraint of fixed number of tetrahedra, we add to the action 



(3.18) a term of the form ejiVa — A'sj, where e is a Lagrange multiplier and is the target number of 
tetrahedra. Essentially, this Lagrange multiplier term modifies the value of the coupling constant fcs. In all 
of the simulations reported below, we have set e — 0.02 and to approximately 10200. These parameter 
values typically yield a one percent variation in over any given ensemble. This number of tetrahedra is 
sufficiently large for physical effects clearly to outweigh finite size effects and sufficiently small for our limited 
computing resources to survey a reasonable portion of the coupling constant space. 

Finally, we tune to a set of coupling constants {fcoi -^i ol, k^[kf), A, a)} on the critical surface of the coupling 
constant space defined by the partition function ( |4.1[ ) for a fixed N^. As our notation suggests, we first select 
values for /cq. A, and a and then tune to the associated critical value k^ of ^3. The critical value k^ is that 
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for which N3 remains approximately constant while a simulation runs. In this sense our model is only well- 
defined at the critical value: for any other value the number of tetrahedra either increases without bound or 
plummets to zero. 

With these conditions established, a simulation begins with the generation of an initial triangulation 
having the topology x composed of N3 tetrahedra. Using the Pachner moves adapted to causal 
dynamical triangulations, as described, for instance, in [S], we run a standard Metropolis algorithm to 
generate an ensemble of spacetimes representative of the weighting defined by the partition function (4.1 1, 
sampling only after a period of thermalization. We sample the representative spacetimes generated every 
one hundred sweeps, a single sweep comprising N3 attempted Pachner moves. Once collected, we estimate 
the expectation values of observables as averages over the ensemble. 

Testing our code is a nontrivial matter: without a known nonperturbative quantization of (2 + 1)- 
dimensional projectable Hofava-Lifshitz gravity, we possess no standard of comparison for our results. Of 
course, this situation also pertains to causal dynamical triangulations formulated with the Regge action. Our 
code is a modified version of that reported in [401 , which has yielded independent corroboration of the results 
of [nmnmilllllllSllIlllHllin]- We have mn the modified code at the general relativistic values of the coupling 
constants A and a to check that we correctly reproduce the results of these references. In figure 3 we present 
depictions of two representative spacetimes — one in the physical phase and one in the decoupled phase of 
causal dynamical triangulations — generated by our code. In figure 4 we plot the ensemble average spectral 
dimensions — discussed further below — for the two ensembles to which the representative spacetimes in figure 
3 belong. Up to finite size effects currently under investigation, these measurements agree quantitatively 
with those of [121 [Ml HDj- In figure 8(a) below we also show the ensemble average discrete 2- volume as 
a function of discrete time for the second of these two ensembles. This data, as well as the fit to it, are 
consistent with the findings of [SlIHl- Beyond these checks, we rely on the plausibility of our new results as 
a test of our code's correctness. 




(a) 



(b) 





10 20 30 40 50 60 



10 20 30 40 50 60 



Figure 3: Depictions of representative spacetimes showing the number A^^^ of spacelike triangles as a function 



of discrete time r. [(a)] Phase A (fco = 6.00, ^3 = 1.85, A = 1.00, a = 0.00) [(b)] Phase C (fco = 1.00, k^. = 0.75, 
A = 1.00, a = 0.00) 



4.2 Phase Diagram 

Our model's coupling constant space is 4-dimensional. Based on the phase structures of both (2 + 1)- and 
(3+l)-dimensional causal dynamical triangulations, we expect that for fixed our model is only well-defined 
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(a) 



(b) 



Figure 4: The ensemble average spectral dimension (dg) as a function of diffusion time a. (a) Phase A 
{ko = 6.00, fca = 1.85, A = 1.00, a = 0.00) [(b)] Phase C {ko = 1.00, fcg = 0.75, A = 1.00, a = 0.00) 



on a 3-dimensional subspace. This critical surface approximates the so-called infinite volume surface, which 
corresponds to the limit in which A'3 increases without bound while the lattice spacing a vanishes such that 
the product N^a^ remains constant. Supposing that our model possesses a second order phase transition, 
at which we could define its continuum limit, this transition must be located at a phase boundary on the 
critical surface. 

Now, a 3-dimensional subspace of largely unknown extent represents a formidably expansive space to 
explore numerically. Accordingly, we have limited our initial investigations to the subspace of the coupling 
constant space consisting of the A — a plane near the origin for fixed fcp. Specifically, we set fcp — 1.00, select 
values for A and a, and then tune to the value of k^ on the critical surface. With A — 1.00 and a = 0.00 this 
corresponds to a point in the physical phase of causal dynamical triangulations for the Regge action. 

Within this region we have generated forty-seven ensembles of representative spacetimes, each for a 
different set {ko, A, a, k^{kQ, A, a)} of the coupling constants. We display in figure 5 the critical surface and 
the associated phase structure as ascertained thus far. Our explorations indicate the existence of three 
phases: a phase contiguous with the physical phase of causal dynamical triangulations for the Regge action 
that we call phase C; a phase emerging for sufficiently large values of A that we call phase D; and a phase 
emerging for sufficiently large values of a that we call phase E. We are not entirely certain that phases 
D and E are distinct: our measurements might instead indicate modulation of a single phase across the 
relevant region in the A — a plane, a possibility under active investigation. On the other hand, we are quite 
certain that phases D and E are not artifacts: our simulations exhibit the characteristic lengthening of the 
autocorrelation time near the phase boundaries, and the geometric properties of phases D and E persist 
under an increase in the total number of spacelike hypersurfaces. The presence of these novel phases of 
extended geometry counts as the first explicit indication that the spatial curvature squared terms in the 
action (3.18) exert a significant effect as opposed to being renormalized to irrelevance. This nicely matches 



the scaling behavior of the curvature squared terms expected from the analytic approach to Hofava-Lifshitz 
gravity. 

In figure 6 we depict spacetimes representative of each of the three phases. Specifically, we plot the 
number A^^^ of spacelike triangles — a measure of the discrete 2-volume — as a function of discrete time r. 
Representative spacetimes in phase C are characterized by a single correlated accumulation of tetrahedra 
spread across a substantial range of r. Representative spacetimes in phase D are characterized by an 
intermittent series of accumulations of tetrahedra each spread across a small range of r. Representative 
spacetimes in the phase E are characterized by a moderately uniform distribution of tetrahedra spread 
across the entire range of r. 



4.3 Physicality and Semiclassicality 

We now present preliminary evidence suggesting that the ensemble average geometries in all three phases 
show signs of being both physical and semiclassical. This evidence stems from analyses of two geometric 
observables associated with our ensembles of representative spacetimes: the spectral dimension as a function 
of diffusion time and the discrete 2-volume as a function of discrete time. Based on these analyses, we draw 
certain comparisons to the relevant solutions of the classical theory now in imaginary time presented in 
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(a) 



(b) 



Figure 5: (a) The critical surface as explored thus far. |(b)| The critical surface projected onto the \-a plane 
showing phases C, D, and E respectively in blue circles, magenta squares, and orange diamonds. 



appendix |Aj 

4.3.1 Spectral Dimension 

The spectral dimension of a space measures its effective dimensionality as experienced by a random walker. 
We determine the ensemble average spectral dimension by the method described, for instance, in jl9l . 
Specifically, we directly measure the return probability Pr{<j) as a function of diffusion time a for each 
representative spacetime in a given ensemble. We then compute the ensemble average spectral dimension as 

(4(.)>^-2^1^ (4.2) 

employing an appropriate discretization of the derivative. 

In figure 7 we display plots of the ensemble average spectral dimension as a function of diffusion time 
for each of the three phases. Before interpreting these plots, we should comment on the evident bifurcation 
for small values of diffusion time. This effect reflects the discrete nature of our triangulated spacetimes: 
random walks of even and of odd lengths yield diverging estimates for the spectral dimension on scales 
sufficiently short in comparison to the discretization scale. Increasing the number of tetrahedra comprising 
each spacetime pushes the bifurcation scale towards smaller diffusion times, revealing the physical nature of 
the spectral dimension on such scales. 

Now compare the plots of figure 7 to those of ordinary causal dynamical triangulations in figure 4. All 
three of the Hofava-Lifshitz spectral dimensions much more closely resemble that of figure 4(b) — for the 
physical phase of causal dynamical triangulations — than that of figure 4(a) — for the decoupled phase of 
causal dynamical triangulations. The resemblance between the plots of figures 7(a) and 4(b) is particularly 
close, as we might expect, since the two ensembles are both in phase C. In particular, the spectral dimensions 
shown in figure 7 all peak at values between 2 and 3 for small diffusion times and then gradually decay for 
large diffusion times. We expect this decay at large diffusion times: these spacetimes are compact and have 
curvature. 

For (2 + l)-dimcnsional causal dynamical triangulations the spectral dimension in the physical phase 
reaches 3 for ensembles characterized by larger values of N-j, HO] . Such measurements provide a key piece 
of evidence for the semiclassical nature of the phase's ensemble average geometry. This suggests that the 
maxima of the spectral dimension plots in figures 4(b) and 7(a) are depressed by finite size effects. Assuming 
that such effects are comparable for our ensembles in phases D and E, we are led to conclude that the 
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(a) 



(b) 
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Figure 6: Depictions of representative spacetimes showing the number A'^f^ of spacehke triangles as a function 
of discrete time r. [(a)] P hase C (fco = 1.00, fcg = 0.77, A = 0.75, a = 0.50) [(b)]Phase D (fco = 1.00, h = 0.73, 
A = 3.50, a = 0.00)](c)1 Phase E (/cq = 1.00, fcg = 0.72, A = 1.00, a = 2.60) 



maxima of the spectral dimension in these phases is slightly below the topological value of 3. In Hofava- 
Lifshitz gravity the spectral dimension is generally predicted to flow under renormalization, so potentially 
the plots of figure 7 reflect such behavior [32] . We are currently running simulations in phases D and E 
for larger values of Ny, in the hope of resolving these issues. Nevertheless, we maintain that the similarities 
in form between the spectral dimensions depicted in flgure 7 and in figure 4(b) hint at the physicality and 
semiclassicality of phases C, D, and E. 

4.3.2 Discrete 2- Volume 

The foliated structure of causal triangulations allows for the measurement of certain quantities as functions 
of the discrete time coordinate. The number iV^^ of spacelike triangles on a given spacelike hypersurface is 
one such quantity. Although the value of N^'^ on a single spacelike hypersurface is not physically meaningful, 
the set {A^^^(St-)}!^^]^ does contain physical information. 

For an ensemble in phase C, there exist straightforward techniques for performing a coherent ensemble 
average of {7V|'-^(S^)}!^^j. (See, for instance, [H].) These techniques rely on the characteristic feature of 
phase C spacetimes: the single accumulation of tetrahedra as depicted in figures 3(b) and 6(a). Intuitively, 
each method functions to align an appropriately defined center of this accumulation with the central value 
of the discrete time coordinate. Once accomplished for all of the representative spacetimes in an ensemble, 
the coherent ensemble average of {A'^-^(St-)}!^^]^ is defined by the discrete timewise average. 

We employ a method that we call equal discrete 2-volume splitting. The name refers to the algorithm 
for appropriately defining the center of the accumulation of tetrahedra. Formally, given a representative 
spacetime, form all divisions Di of that spacetime into two sets each of spacelike hypersurfaces, maintaining 
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(a) 



(b) 



(c) 



Figure 7: The ensemble average spectral dimension {dg) as a function of diffusion time a. (a) Phase C 
{ko = 1.00, /c3 = 0.77, A = 0.75, a = 0.50) [(b)] Phase D (/cq = 1-00, h = 0.73, A = 3.50, a = 0.00) |(c)| Phase 
E {ko = 1.00, ks = 0.72, A = 1.00, a = 2.60) 



the periodic discrete time order of the spacelike hypersurfaces: 



D2 = {{s2,...,St+i},{e.+2,...,Si}} 



(4.3a) 
(4.3b) 



Dt = {{St,...,St_i},{Et,...,St-i}} 



(4.3c) 



Then select the two particular divisions D^q and D^q that most nearly equalize the discrete 2-volume summed 
over the spacelike hypersurfaces in each set of the division: the quantity 



E 



E 



(4.4) 



1} 



is minimized for both D(.q and -Deg- 



Note that De„ and D^q only differ in the order of their two sets. Next 



relabel the discrete time coordinate over the set of values {- 



T 
2 ' 



+ 1,.. 



T 
2 



1, 4} so that the values 



{— 1^, . . . , — i} label the first set of spacelike hypersurfaces and the values {i, . . . , ^} label the second set of 
spacelike hypersurfaces in both D^q and Deq- Finally choose the division for which the values of A^^^(I]_i) 
and N2^{T,i) are greatest. With all of the spacetimes in an ensemble so aligned, perform a discrete timewise 

average of iV^^ over the representative spacetimes. 

In figure 8 we show the results of the equal discrete 2-volume splitting average for two ensembles in phase 
C. The data points indicate the equal discrete 2-volume splitting average value of {Ni^{^r)}r=i; the light 
vertical bars indicate one standard deviation of error. The thin curve is a one parameter fit of the data 
points within the central accumulation of tetrahedra to the functional form 



ir) 



■ cos 



(4.5) 



which is a discretization of the 2-volume as a function of global time of Euclidean de Sitter spacetime. (See 
[24] and [H] for a derivation of this discretization in 2 -I- 1 and 3 4- 1 dimensions, respectively.) Substantial 
evidence already demonstrates that, for the general relativistic values of A and a, the ensemble average 
geometry in phase C closely approximates Euclidean de Sitter spacetime [H [51 [TOl [HI [HI [13 [HI HH [IHl [21] ■ 
The plot in figure 8(b) provides the first evidence that the ensemble average geometry also has this property 
for the broader range of nonrelativistic values of both A and a. 
Hofava-Lifshitz gravity admits Euclidean de Sitter spacetimes 



for A > ^ and a = 0, the solutions (A.IO) 



and ( A.ll ) coincide in their description of such spacetimes. Note, however, that for relatively small values of 
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Figure 8: The ensemble average number (A'^^) of spacelike triangles as a function of discrete time r. (a) 
Phase C (fco = 1.00, h = 0.75, A = 1.00, a = 0.00) with So = 0.46 [(b)] Phase C (fco = 1.00, h = 0.77, 
A = 0.75, a = 0.50) with Sq = 0.48 



Figure 9: The ensemble average number (iV^^) of spacelike triangles as a function of discrete time t for four 
ensembles in phase C differing in their respective values of the coupling constant A: blue circles (fco = 1.00, 
fca = 0.85, A = -1.00, a = 0.00), magenta squares (fco = 1-00, fcs = 0.79, A = 0.00, a = 0.00), orange 
diamonds (fco = 1.00, fcg = 0.75, A = 1.00, a = 0.00), purple triangles (fco = 1.00, h = 0.74, A = 2.00, 
a = 0.00). 



a, the solutions (A. 10 1 and (A. 11 1 do not deviate too markedly from Euclidean de Sitter spacetime. Given 
the presence of finite size effects and the inherent error in our measurements, one of these solutions for a ^ 
may well fit better the ensemble average geometry in phase C. Our current data are not sufficiently detailed 
for us to make a conclusive statement. Furthermore, there exists some evidence [19] that the ensemble 
average geometry in (2 + l)-dimensional causal dynamical triangulations is a deformed version of Euclidean 
de Sitter spacetime. Supposing that this is the case, the process of taking the continuum limit must generate 
some additional terms besides those in the Regge action. 

The two plots in figure 8 quite closely resemble one another — in both spatial and temporal extents of the 
central accumulation of discrete 2-volume — even though that of figure 8(a) is for general relativistic values 
of A and a while that of figure 8(b) is for nonrelativistic values of A and a. To dispel the suspicion that 
the and i?2 terms are renormalized to irrelevance in phase C, we present in figure 9 four plots of the 
equal discrete 2-volume splitting average of {N2^{T,t-)}'^^i each for a different value of A. Clearly, varying 
A affects the ensemble average geometry that emerges. 

For phases D and E we currently do not know how to coherently average {N2^{'Sr)}T=i over an ensemble. 
In an effort to determine a method, we computed the ensemble average power in the discrete Fourier transform 
of {N2^{T,T)}'^^-^^, the results of which we display in figure 10. As these plots show, there is no notable 
periodicity present in these ensemble's average geometry since virtually all of the power falls in the zero 
frequency mode. This lack of periodicity may in fact be indicative of the semiclassical nature of phase E. 
With a > the classical equations of motion (A.l) and (A. 2 1 are also satisfied for the Ansatz (A. 4) when 
the squared scale factor a^(t) has the constant value The depiction in figure 6(c) of a representative 

spacetime in phase E resembles to a certain extent a spacetime with approximately constant scale factor, 
and the plot of figure 10(b) reinforces this interpretation. Moreover, this accords with our expectation that 
as currently devised the Markov chain Monte Carlo simulations converge on spacetimes globally minimizing 
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the action (3.18). 
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Figure 10: The ensemble average power (|cp) as a function of discrete frequency v in the discrete Fourier 
transform of the number -/V^^ of spacelike triangles as a function of discrete time r. (a) Phase D {ko = 1.00, 
ks = 0.73, A = 3.50, a = 0.00) [(b)] Phase E {ko = 1.00, k^ = 0.72, A = 1.00, a = 2.60) 



At first glance, however, the deviations from constant scale factor as illustrated in figure 6(c) do not appear 
particularly small. To make a proper judgment, we compare the deviations in {7V^^(I]T-)}i^^]^ between the 
ensemble in phase C associated with figure 6(a) and the ensemble in phase E associated with figure 6(c). 
Specifically, for each representative spacetime of this ensemble in phase E, we calculate the discrete time 
averaged deviation A^si, of the discrete 2-volume from the mean. The ensemble average (A^si) provides a 
reasonable measure of the deviations from a time-independent geometry. For this ensemble in phase C, we 
calculate the deviations (A^si,)mm and (A^sj:,)^^^, of the discrete 2-volume from the mean for the spacelike 
hypersurfaces of minimal and of maximal {N2^) within the central accumulation of tetrahedra. The values 
(A^si)„ij„ and {Aj^sL)„iax provide a reasonable measure of the range of deviations from a Euclidean de 
Sitter geometry. We find that 



\^^2 /max 



0.20 



< 



0.69 



(A 



< 



and that 



(A 



= 0.35 > 



= 0.27 > 



\-'*2 Irmn 



0.78 



c 



3/min 



= 0.22. 



(4.6) 



(4.7) 



To achieve a more faithful comparison, we have considered the above two ratios instead of the deviations 
themselves. Although not definitive, the fact that the ratios for phase E fall between those for phase C lends 
credence to our interpretation of the semiclassical nature of phase E. 

Supposing that further analysis supports the conjecture that time-independent geometries dominate 
phase E, how could one then interpret the C-E phase transition? A possible answer involves the global 



gravitational Hamiltonian 'H± of (2.101. Recall that, in our implementation of Hofava-Lifshitz gravity into 



the framework of causal dynamical triangulations, we chose not to impose the condition (2.9 1 as a constraint. 



instead treating 'H± as a conserved global charge, essentially the total energy. In our statistical ensemble 
whether Hj^ vanishes or has a nontrivial spectrum on physical states thus becomes a question of dynamics. 
If the dynamics enforce the vanishing of H_l on physical states, then we might interpret this phenomenon 
as dynamical confinement of the gravitational charge. A similar phenomenon has recently been discovered 
in (2 -f l)-dimensional relativistic chiral gravity [44]. There not just one gravitational charge but an infinite 
hierarchy of conserved chiral charges are confined, that is, dynamically vanish on all physical states of finite 
energy. 

Now consider the hypothetical phase in which is dynamically confined. On spacetime geometries of the 
Friedmann-Lemaitre-Robertson- Walker type, the Hamiltonian constraint equation becomes the Friedmann 
equation for the scale factor. The Friedmann equation precludes the ground state geometry from being 
time-independent, forcing a cosmological, de Sitter-like evolution of the universe. In the context of causal 
dynamical triangulations, we observe this behavior in phase C. 
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In the hypothetical deconfined phase, on the other hand, the situation is more reminiscent of a typical 
condensed matter system: the total Hamiltonian simply measures the system's energy levels with the ground 
state identified as the lowest energy state, typically static. We apparently observe this behaviour in our 
phase E. Accordingly, we speculate that the C-E phase transition may be viewed as the deconfinement of 
the gravitational charge H±. 

5 Outlook 

Motivated by a suite of striking similarities between Hof ava-Lifshitz gravity and causal dynamical triangula- 
tions, we constructed a discrete version of the action for (2 + l)-dimensional Hofava-Lifshitz gravity adapted 
to the formalism of causal dynamical triangulations. Using this action in Markov chain Monte Carlo simula- 
tions of the corresponding path integral, we found significant evidence for the existence of extended phases 
of geometry approximating certain classical solutions of Hofava-Lifshitz gravity. Quantum Hofava-Lifshitz 
gravity thus appears amenable to and compatible with the techniques of causal dynamical triangulations. 

Since we have only just initiated this study, the prospects for further research are expansive. In the 
near term we have three primary goals for ongoing research. First, we plan to map more extensively our 
model's coupling constant space to determine the extent of the phases so far discovered and to identify 
the locations of phases not yet discovered. Specifically, we hope to ascertain how the curvature squared 
terms affect the decoupled phase of causal dynamical triangulations. Second, we wish to further our analysis 
of the semiclassical natures of phases C, D, and E along the lines of [H UHl Uni HI]- Third, we want to 
establish the orders of our model's phase transitions. In this direction we have attempted to determine order 
parameters for the C-D and C-E phase transitions; unfortunately, none of our trial order parameters have 
yet distinguished between these adjacent phases. If our speculation about the confinement-deconfinement 
nature of the C-E phase transition is correct, then the appropriately defined ground state energy may act as 
an order parameter. 

In the long term the richness of the literature on Hofava-Lifshitz gravity provides a host of directions 
for continuing research. First of all, recall that the simplest versions of (2 -I- l)-dimensional Hofava-Lifshitz 
gravity, both projectable and nonprojectable, possess a single propagating scalar degree of freedom. This is 
in stark contrast to the case of not only (2 + l)-dimensional general relativity, but also the nonrelativistic 
generally covariant version of Hofava-Lifshitz gravity constructed in |33j . This scalar mode has generated 
considerable controversy with regard to its phenomenological viability [461 152] , but its full dynamics, espe- 
cially around the stable ground state, remains poorly understood. The (2 + l)-dimensional theory should 
provide a clear window into the scalar mode's dynamics since the complications of propagating tensor modes 
are absent. By studying the nonperturbative dynamics of the scalar mode using causal dynamical trian- 
gulations, we hope to shed light on this issue. Our first challenge is the identification of an observable in 
causal dynamical triangulations that measures the number of local propagating degrees of freedom or that 
at least distinguishes between the absence and presence of local propagating degrees of freedom. Given the 
relative ease of extracting correlation functions of observables associated with the spacelike hypersurfaces of 
causal triangulations, we are currently working to understand the behavior of our model's conformal mode 
in relation to the scalar mode. We hope that such an investigation might illuminate the nature of the scalar 
dynamics. 

Relatedly, perturbations of projectable Hofava-Lifshitz gravity about (3 + l)-dimensional Minkowski 
spacetime generate instabilities |41| . This finding is of course not surprising: Minkowski spacetime is not 
the ground state of the commonly analyzed models. The situation is much improved on the background of 
de Sitter spacetime [25l[53]. As the coupling constant A approaches unity, however, the higher derivative 
terms become relevant, leading to a breakdown of the linearized analysis. The authors of [45l [53] have 
suggested that a nonperturbative Vainshtein mechanism might take effect, rendering this limit continuous 
to the general relativistic values of the coupling constants. Potentially, we could uncover this transition 
behavior in Markov chain Monte Carlo simulations. Indeed, the critical surface depicted in figure 5(a) 
provides a small but intriguing piece of evidence for that possibility: the apparent geometric feature along 
the A = 1 direction. 

To assess the phenomenological viability of at least the simple version of Hofava-Lifshitz gravity on which 
we have focused our study, we need to study the renormalization group flow of the coupling constants for the 
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purpose of comparing the long distance behavior to that of general relativity. In [35] Henson proposes a coarse 
graining procedure for causal dynamical triangulations on which one could try to base a renormalization 
group procedure. If this scheme proves apt - a question that we are currently exploring - then we plan to 
employ the procedure to study the renormalization group flows of our model. Ultimately, we would like to 
explore the general conjecture formulated by several groups PI IM] suggesting that Hofava-Lifshitz gravity 
and causal dynamical triangulations belong to the same universality class. This circumstance would neatly 
explain the remarkable resemblances between these two approaches to the quantization of gravity. 
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A On the Classical Equations of Motion 



We collect here the classical equations of motion stemming from the action (2.81 for constant lapse and their 



solutions relevant to the Markov chain Monte Carlo simulations reported above. In Riemannian signature 
and for spacelike hypersurfaces having the topology of 5^, the equations of motion are 



= V,^''J(i,x) 



(A.l) 



and 




(i , x) J + -Y' {t, x) [Kk, {t, x)X^' (t , x) - XK^ (t , x) - aRj (i , x) + 2A] 



where 



-2if''(t,x)X/(i,x) + 2Aif(t,x)/s:*J'(t,x) +2aV''V^i?2(t,x) - 2aV^i?2(t, x)7^^'(t, x) 

+ i [\/iN\t, x)^-'"' {t, x) + ViiV^' (<, x)^^'(i, x) - tt'^ {t, x) ViiV'(t, x)] (A.2) 

Tr,j{t,x) = if„-(t,x) - AX(t,x)7,,(t,x). (A.3) 



(A.2) result from variation of the spatial metric tensor. Again we exclude the nonlocal integral constraint 



The momentum constraint ( A.l ) results from variation of the shift vector, and the metric equations of motion 



(|2.9 ) arising from variation of the lapse. 

We seek solutions to the equations of motion ( |A.1[ ) and ( |A.2p in the form of the Friedmann-Lemaitre- 



Robertson- Walker Ansatz 



(A.4) 



for spatially homogeneous and isotropic spatial metric tensor "f{t, x) and identically vanishing shift vector 
N(i,x). Here, 7(x) is the metric tensor on the round 2-sphere, and a{t) is the scale factor. Under these 



assumptions the momentum constraint (A.l I is trivially satisfied. The metric equations of motion (A.2) are 
satisfied if and only if 



a 



1 



A 



which implies that 



d_ 
dt 



(2A- 1)C 
2A 



1 d^ajt) 

^ dt2 ~ 2(2A - 1) a*{t) ^ 2A- 1 

C2 (2A - 1) 



2a 
2A - 1 



A 



4A 
2A- 1 



a'it)- 



C(2A- 1) 
2A 



(A.5) 



(A.6) 
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for constant of integration C. Letting 



(A. 6 1 becomes 



which has the solution 



u{t) = a^{t) 



C(2A-1) 
2A 



du{t) 
dt 



4A 
2A- 1 



u\t) 



2a 



2A- 1 



(2A- 1) 
A ' 



uit) 



ic^ {2x-iy 

4A2 



2A 



cos 2 



A 

2A- 1 



t + S 



for a second constant of integration S. In terms of a{t), the solution (A. 9 1 is 

C(2A-1) /C2 (2A-1)^ 



<W = 



2A 



4A2 



2A 



cos 2 



A 



2A- 1 



t + S 



(A.7) 



(A.8) 



(A.9) 



(A.IO) 



for positive values of a and 



2A 



C(2A-1) /C2(2A-1)^ |a| 



4A2 



2A 



cos 2 



A 



2A- 1 



t + S 



(A.ll) 



for negative values of a. We assume that A > as in the previous studies of causal dynamical triangulations. 
Note that the solution (A.ll) has a finite time extent dictated by the zeros of a^{t) whereas the solution 
(A.lOl has no such restrictions. Also note that for vanishing a the solutions (A.IO) and (A.ll) both reduce 
to that of Euclidean de Sitter spacetime. 



B On the Geometry of Causal Dynamical Triangulations 

For the tetrahedra depicted in figure 1, the Lorentzian dihedral angles between spacelike and tinielike faces 
are 



(3,1) 


TT 


i\og 


(2,2) 


= i log 




(1,3) 


TT 

- 2^ 


i\og 



and the Lorentzian 3-volumes are 



l + 2v/3^7TT 
3-2^/2V2^ 

4:7]+ I 



1 + 2^3^7+1 



(B.la) 
(B.lb) 
(B.lc) 



V, 



t'^ = ^V^V+la' (B.2a) 



V, 



(1.3) _ 



L 



^V3^+la3. (B.2c) 
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As previously mentioned, Wick rotation consists in analytically continuing r] in the lower half complex 
plane [8 . If the argument of a square root becomes negative as a result of the Wick rotation, then we replace 
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it by the negative of the argument multiphcd by ~i. The Lorentzian dihedral angles (B.l ) arc thus continued 
to their respective Euclidean values 



q(3,l) 



q(2,2) 
^E 



q(l,3) 
^E 



- — cos ' 



-1 / 4ry-3 



— — cos ' 



and the Lorentzian S-volumes are thus continued to their respective Euclidean values 



(B.3a) 
(B.3b) 
(B.3c) 



V, 



(2,2) 



12 

i 



6V2 



V, 



(1,3) 



12 



(B.4a) 
(B.4b) 
(B.4c) 



Note that the magnitude of rj must be greater than ^ to prevent the tetrahedra from becoming degenerate. 
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